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^  The  infljence'of  the  ocean-atmosphere  interface  on  the  radiance  distribu¬ 
tion  in  both  the  ocean  and  atmosphere  is  investigated.  At  visible  wavelengths 
in  the  real  ocean  just  below  the  surface,  the  downwelling  radiance  is  a  maximum 
within  the  critical  angle  and  drops  off  by  one  or  two  orders  of  magnitude  to¬ 
ward  the  horizon.  The  usual  explanation  that  this  is  due  to  the  sky  radiation 
concentrated  within  the  critical  angle  and  the  total  internal  reflection  of  the 
weak  upward  radiance  at  the  ocean  surface  at  angles  outside  the  critical  angle 
is  too  simplistic.  There  are  two  other  important  factors:  (0  the  atmosphere^ 
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20.  Abstract 


must  have  sufficient  optical  thickness  so  that  appreciable  sky  (multiple 
scattered)  radiation  develops  at  all  angles  in  the  downward  direction;/^  D 
C2j  the  water  must  have  appreciable  absorption  so  that  the  upwelling 
radiance  just  below  the  ocean  surface  is  only  a  small  fraction  of  the 
downwelling  radiance  entering  the  ocean ~  Examples  show  that,  if  either 
one  of  these  conditions  is  not  satisfied,!  an  entirely  different  radiance 
distribution  develops.  The  variation  of  Ahe  following  quantities  with 
depth  is  studied:  radiance,  vector  and  scalar  irradiance,  distribution 
function,  reflectance,  and  heating  rateJ^The  radiance  distribution  in 
a  homogeneous  medium  is  compared  with  that  in  the  same  medium  with  an 
atmosphere-ocean  interface  at  various  depths  within  the  medium.  Most 
of  the  calculations  are  done  for  Rayleigh  scattering  centers,  but  some 
results  are  given  for  Mie  type  scattering. 
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ABSTRACT 

The  influence  of  the  ocean-atmosphere  interface  on  the  radiance 
distribution  in  both  the  ocean  and  atmosphere  is  investigated.  At 
visible  wavelengths  in  the  real  ocean  just  below  the  surface,  the  down- 
welling  radiance  is  a  maximum  within  the  critical  angle  and  drops  off 
by  one  or  two  orders  of  magnitude  toward  the  horizon.  The  usual  ex¬ 
planation  that  this  is  due  to  the  sky  radiation  concentrated  within 
the  critical  angle  and  the  total  internal  reflection  of  the  weak  upward 
radiance  at  the  ocean  surface  at  angles  outside  the  critical  angle  is 
too  simplistic.  There  are  two  other  important  factors:  1.  the  atmo¬ 
sphere  must  have  sufficient  optical  thickness  so  that  appreciable  sky 
(multiple  scattered)  radiation  develops  at  all  angles  in  the  downward 
direction;  2.  the  water  must  have  appreciable  absorption  so  that  the 
upwelling  radiance  just  below  the  ocean  surface  is  only  a  small  fraction 
of  the  downwelling  radiance  entering  the  ocean.  Examples  show  that.  If 

either  one  of  these  conditions  is  not  satisfied,  an  entirely  different 
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radiance  distribution  develops.  The  variation  of  the  following  quqHtfV  -- 
ties  with  depth  is  studied:  radiance,  vector  and  scalar  irradianqe,  ... . 
distribution  function,  reflectance,  and  heating  rate.  The  radiance 
distribution  in  a  homogeneous  medium  is  compared  with  that  in  the  s*me 
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medium  with  an  atmosphere-ocean  interface  at  various  depths  within  the., 
medium.  Most  of  the  calculations  are  done  for  Rayleigh  scattering 
centers,  but  some  results  are  given  for  Mie  type  scattering. 
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I.  Introduction 

Divers  marvel  at  the  beautiful  patterns  of  light  visible  within 
the  ocean.  When  the  surface  is  rough,  it  is  difficult  for  the  eye  to 
provide  information  for  a  time-averaged  radiance  from  the  dancing 
points  of  light.  On  the  other  hand  when  the  surface  is  calm,  the  entire 
world  above  the  ocean  is  condensed  by  refraction  of  the  rays  at  the 
surface  into  a  cone  extending  from  the  zenith  to  the  critical  angle. 

If  the  surface  is  perfectly  smooth,  the  radiance  drops  abruptly  by  one 
or  two  orders  of  magnitude  as  the  viewing  angle  moves  outside  the 
critical  angle.  The  usual  explanation  is  that  the  sky  radiation  is 
concentrated  within  the  critical  angle  in  the  ocean,  while  the  down- 
welling  radiation  outside  the  critical  angle  is  derived  from  the  total 
internal  reflection  of  the  much  weaker  upwelling  radiation  just  beneath 
the  ocean  surface. 

Although  this  explanation  is  correct  as  far  as  it  goes,  it  does 
not  point  out  the  dependence  on  several  other  properties  of  the  medium. 
The  influence  of  the  absorption  of  the  water  in  reducing  the  upwelling 
radiance  within  the  ocean  is  crucial,  as  is  the  existence  of  a  suffi¬ 
ciently  deep  atmosphere  so  that  appreciable  downwelling  sky  radiance 
develops.  The  dependence  of  the  radiance  distribution  on  these  and  the 
other  optical  properties  of  the  medium  is  thoroughly  explored  in  the 
following  sections.  For  ease  of  calculation  all  results  given  In  Sect. 
II  -  VII  are  for  Rayleigh  scattering.  The  effect  of  M1e»type  scattering 
is  Illustrated  in  Sect.  VIII.  A  smooth  ocean  surface  Is  assumed  for 
all  these  calculations. 
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II.  Method  of  Calculation 

The  Matrix  Operator  Method  was  used  for  all  calculations  reported 

here.  Theoretical  expositions  and  applications  of  the  theory  to  an 
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inhomogenous  medium  have  been  previously  presented  by  us.  The 
study  of  radiative  transfer  problems  pertaining  to  the  atmosphere-ocean 
system  requires  that  the  theory  be  modified  to  include  reflection  and 
refraction  at  the  ocean's  surface. 

A  new  fundamental  problem  was  encountered  as  a  result  of  the  attempt 
to  implement  these  modifications.  The  problem  encountered  is  a  direct 
result  of  the  law  of  refraction.  In  the  standard  matrix  operator  for¬ 
malism,  the  set  of  zenith  angles  used  for  the  purpose  of  calculation 
is  derived  from  a  quadrature;  usually  a  set  of  abcissas  equal  to  the 
positive  points  of  either  a  Gauss  or  Lobatto  quadrature.  Within  the 
modified  theory,  once  a  set  of  zenith  angles  is  selected  to  represent 
either  the  atmosphere  or  the  ocean,  zenith  angles  (quadrature  points) 
for  the  other  region  are  specified  by  Snell's  Law.  No  proven  way  of 
;  assigning  weights  to  the  non-linearly  mapped  points  was  available  to 

>,  >  us.  For  this  reason,  selection  of  the  atmospheric  quadrature,  selection 

%:  of  the  ocean  quadrature  and  formulation  of  the  mapping  relating  the 

y<  _"Av- 

m',-,  weights  of  these  quadratures  was  crucial  to  the  Implementation  of  the 

*  /  J  * 

modified  theory.  All  other  problems  related  to  inclusion  of  the  di¬ 
electric  interface,  i.e.  construction  of  reflection  and  transmission 
operators  for  the  interface  and  phase  function  normalization  criteria 
were  attendant  to  or  subordinant  to  this  problem. 


Four  constraints  arose  naturally  In  the  process  of  selecting  the 
quadratures  to  represent  the  atmosphere-ocean  system.  First  the  quad¬ 
rature  for  the  atmosphere  had  to  normalize  a  phase  function  adequate 
to  represent  scattering  by  atmospheric  aerosols.  Second,  the  ocean 
quadrature  had  to  normalize  the  highly  assymetrlc  phase  functions  used 
to  represent  hydrosol  scattering.  Third,  the  related  quadratures  had 
to  guarantee  conservation  of  energy  transmitted  and  reflected  by  the 
dielectric  Interface.  Finally,  both  the  selected  quadratures  had  to 
retain  sufficient  strength  to  Integrate  properly  the  intensity  distri¬ 
butions  arising  In  their  respective  regions. 

Our  solution  to  this  problem  was  to  use  a  set  of  M  mapped  Gauss 
points  p .  as  our  atmospheric  quadrature.  The  quadrature  points  5. 

J  J 

within  the  ocean  were  then  found  by  a  simple  application  of  Snell's 
Law,  namely 


“j  *  ['  -  c 


in 


or 


Cj  *  [1  *  (1  *  Uj2) n‘2]v2  ,  (2) 

where  pj 
face,  Cj 
n  is  the 


are  the  cosines  of  the  angles  of  incidence  from  above  the  inter 
are  the  cosines  of  the  refracted  angles  within  the  ocean,  and 
index  of  refraction  of  water  relative  to  air.  The  weights  for 


these  ocean  quadrature  points  within  the  acceptance  cone  are  found  by 


(3) 


CJ  *  <V«J"  )CJ» 

where  C*?  is  the  calculated  ocean  quadrature  weight  corresponding  to  t;. 

J  J 

and  is  the  mapped  Gauss  weight  of  the  corresponding  atmospheric 
quadrature  point  p..  An  additional  set  of  K  Gauss  points  was  mapped 

J 

into  the  regime  from  the  critical  angle  to  the  horizon.  A  complete 

explanation  of  the  reasoning  behind  these  selections  was  presented 
4 

by  Kattawar,  et  al_. 

With  the  quadrature  types  selected,  solutions  to  all  other  problems 
were  then  simple  to  formulate.  The  minimum  order  of  the  atmospheric 
quadrature  was  determined  by  requiring  it  to  allow  the  normalization 
of  a  chosen  aerosol  phase  function  to  be  within  0.05%.  At  the  same 
time,  the  related  ocean  quadrature  was  required  to  normalize  the 
necessary  hydrosol  phase  function  to  similiar  accuracy.  To  minimize 
the  size  of  the  ocean  quadrature  necessary  to  normalize  the  hydrosol 
phase  function,  the  delta  function  approximation  used  by  Potter  was 
employed.5 

With  suitable  quadratures  selected,  the  problem  of  constructing 
reflection  and  transmission  operators  to  represent  the  interface  pro¬ 
perties  became  straightforward.  For  a  scalar  calculation  they  are 
merely  diagonal  matrices  formed  from  the  fresnel  reflection  and  trans¬ 
mission  coefficients,  r(p.)  and  t(p.)f  for  unpolarized  radiation*  and 

J  J 

a  geometric  factor  derived  from  the  change  in  solid  angle  on  entry 
to  or  exit  from  the  ocean.5  The  reflection  and  transmission  operators 
applicable  to  light  traversing  the  interface  from  the  atmosphere  into 
the  ocean  may  be  written: 


(4) 


[R(A,W)]jk  -  RtA.Wjj  5jk  , 

[T(A,W)]jk  =  n2T(A,W)j  6jk  (5) 

for  j  =  1  to  M. 

For  light  traversing  the  interface  from  the  ocean  into  the  atmosphere 
the  reflection  and  transmission  operators  can  be  written: 


[R(W,A)].k  =  R(W,A)fijk  , 

(6) 

[T(W,A)]jk  =  rT2T(W,A)..  6jk  , 

(7) 

for  j  =  1  to  M, 


[R(W,A)]j'k'  =  6j 'k '  * 

(8) 

[T(W,A)]r|r  =  0  , 

(9) 

for  j'  *  M  +  1  to  M  +  k. 

Completely  new  computer  algorithms  were  implemented  to  accomodate 
the  modifications  to  the  method.  The  new  algorithms  perform  both  the 
co-decomposition  of  aerosol  and  hydrosol  phase  functions,  and  the  radia¬ 
tion  field  calculations  for  a  combined  inhomogeneous  atmosphere-ocean 
system.  In  addition  facilities  to  allow  calculation  for  any  ocean 
bottom  with  a  known  reflectance  have  been  built  into  the  program.  Two 
special  features  of  particular  value  when  studying  radiative  transfer 


problems  associated  with  the  atmosphere -ocean  system  were  also  built  into 
the  program.  The  first  allows  the  spectra  from  a  complete  wavelength  scan 
to  be  calculated  in  a  single  run.  The  second  allows  calculations  for 
multiple  geometries  to  be  made  within  a  single  run.  The  geometries  of 
the  second  feature  are  limited  to  those  which  can  be  constructed  by  re¬ 
positioning  the  ground  at  the  levels  of  any  detectors  imbedded  in  the 
ocean.  Each  aspect  of  the  new  code  was  tested  by  comparison  to  previous 
calculations^’7”^. 

III.  Effect  of  Water  Boundary 

The  radiance  underwater  has  a  distinct  pattern  well  known  to  under¬ 
water  swimmers  and  divers.  When  looking  upward  at  the  water  surface, 
there  is  a  bright  pattern  of  light  visible  at  zenith  angles  up  to  the 
critical  angle,  approximately  48°.  By  comparison  the  ocean  appears  very 
dark  when  viewed  at  zenith  angles  greater  than  the  critical  angle.  This 
distribution  is  observed  when  the  water  surface  is  calm  (waves  smooth 
out  the  sharp  drop  in  the  radiance  at  the  critical  angle^),  the  bottom 
is  not  too  near  the  surface,  and  the  water  is  reasonably  clear. 

This  usual  underwater  radiance  distribution  has  been  explained 
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many  times  in  the  literature  ’.  All  of  the  light  from  the  sun  and 
sky  that  enters  the  ocean  is  refracted  into  zenith  angles  less  than  the 
critical  angle,  while  no  light  from  the  atmosphere  can  enter  a  calm 
ocean  at  greater  angles.  The  upwelling  light  in  the  ocean  just  below 
the  surface  is  always  small  in  comparison  to  the  downwelling.  If  this 
weak  upwelling  light  encounters  the  ocean  surface  at  angles  greater  than 


the  critical  angle  it  is  totally  internally  reflected  and  thus  becomes 
downwelling  light.  This  is  the  only  source  of  downwelling  light  at 
angles  greater  than  the  critical  angle  and  thus  the  downwelling  light 
is  weak  outside  the  cone  defined  by  the  critical  angle. 

Although  this  explanation  is  correct  as  far  as  it  goes,  some  de¬ 
tails  concerned  us.  A  fundamental  investigation  of  the  manner  in  which 
the  light  field  develops  underwater  seemed  of  interest.  It  soon  became 
obvious  that  the  situation  is  indeed  more  complicated  in  that  other 
factors  also  control  the  light  distribution.  This  is  illustrated  in 
Fig.  1,  which  shows  the  downward  diffuse  radiance  just  below  the  ocean 
surface  as  a  function  of  the  zenith  angle  of  observation  for  a  solar 
zenith  angle,  eQ  =  11.44°.  All  radiance  curves  given  in  this  article 
have  been  averaged  over  all  azimuthal  angles.  There  is  no  absorption 
in  the  atmosphere  unless  otherwise  stated.  Only  Rayleigh  scattering 
is  assumed  in  both  the  atmosphere  and  ocean  unless  mentioned  otherwise. 
This  assumption  greatly  simplified  the  calculations  and  makes  it  easier 
to  study  the  other  factors  that  determine  the  underwater  radiance  dis¬ 
tribution.  The  incoming  solar  radiation  is  normalized  to  unit  flux 
through  a  surface  perpendicular  to  the  beam. 

The  top  curve  in  Fig.  1  is  for  a  set  of  parameters  that  corresponds 
to  the  actual  ocean-atmosphere  system  at  most  visual  wavelengths:  the 
single  scattering  albedo  (fraction  of  radiation  actually  scattered  for 


each  photon  collision)  of  the  ocean,  u>Q  =  0.5  and  the  optical  depth  of 
the  atmosphere  is  taken  as  t  =  0.5.  The  ocean  is  assumed  to  be  ef¬ 


fectively 


of  infinite  depth.  This  radiance  curve  has  the  typical 


variation  observed  in  the  visible  in  the  real  ocean:  the  radiance  is 
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relatively  large  for  zenith  angles  of  observation  less  than  the  critical 
angle  and  small  at  larger  angles.  There  is  a  rapid  change  in  the  radiance 
near  the  critical  angle. 

If  we  merely  change  the  single  scattering  albedo  of  the  ocean  from 
“0  =  0.5  to  wo  =  1  and  leave  all  the  other  parameters  the  same,  an  en¬ 
tirely  different  radiance  curve  is  obtained,  the  bottom  curve  of  this 
figure.  With  this  change  in  the  ocean  albedo,  the  upwelling  radiation 
is  so  much  larger  that,  after  total  internal  reflection  at  the  water 
surface,  the  downwelling  radiation  observed  beyond  the  critical  angle 
is  now  much  larger  than  the  downwelling  radiation  from  the  sun  and  sky 
observed  within  the  critical  angle.  Thus,  the  observed  radiance  curve 
in  the  ocean  depends  on  the  fact  that  the  ocean  water  is  a  fairly  good 
absorber  of  radiation  at  all  visible  wavelengths.  If  it  were  not,  the 
radiance  would  have  a  maximum  at  angles  greater  than  the  critical  angle. 

Let  us  change  just  one  other  parameter  in  our  original  calculation. 
Instead  of  having  a  fairly  substantial  atmosphere  above  the  ocean  sur¬ 
face,  let  us  have  only  a  very  thin  atmosphere,  Tatm  =  0.0001.  The 
radiance  curve  in  the  middle  of  the  figure  is  obtained  for  this  situa¬ 
tion.  Again  the  largest  radiance  values  occur  at  angles  of  observation 
beyond  the  critical  angle.  These  curves  show  only  the  diffuse  radiance, 
i.e.  photons  which  have  been  scattered  one  or  more  times  from  the  direct 
solar  beam.  In  this  case,  there  are  only  a  few  photon  collisions  in  the 
very  thin  atmospheric  layer  and  there  is  virtually  no  sky  radiation  en¬ 
tering  the  ocean.  Thus,  the  downwelling  light  observed  just  below  the 
ocean  surface  is  largely  light  from  the  upwelling  radiation  that  has 
been  internally  reflected  at  the  water  surface  and  thus  changed  into 
downwelling  radiation.  At  angles  greater  than  the  critical  angle,  all 


of  this  radiation  is  reflected  by  total  internal  reflection.  At  angles 
less  than  the  critical  angle,  the  fraction  reflected  can  be  easily  de¬ 
termined  from  the  usual  laws  for  reflection  and  is  found  to  decrease 
rapidly  as  the  angle  of  observation  decreases  from  the  critical  angle. 

In  this  case,  even  though  the  ocean  still  has  -  0.5,  the  radiance 
distribution  is  completely  opposite  from  the  observed  one. 

The  conclusion  from  this  brief  example  is  that  the  observed  radiance 
distribution  in  the  ocean  requires  both  (1)  an  atmosphere  of  reasonable 
optical  thickness  in  order  to  develop  appreciable  sky  radiation  as  well 
as  (2)  an  ocean  that  absorbs  an  appreciable  fraction  of  the  photons. 

If  either  one  of  these  requirements  is  not  satisfied,  the  radiance  dis¬ 
tribution  may  be  the  complete  opposite  of  the  one  usually  observed 
in  the  ocean. 

IV.  Dependence  of  Downwelling  Radiance  on  Ocean  Albedo  and  Atmospheric 

Thickness 

More  detailed  results  are  given  in  this  section  in  order  to  show 
the  dependence  of  the  radiation  field  on  the  ocean  albedo  and  the  at¬ 
mospheric  thickness.  For  simplicity  the  Rayleigh  phase  function  was 
used  for  all  scattering  events  in  both  the  atmosphere  and  ocean  in  this 
section.  In  all  cases  the  ocean-atmosphere  interface  reflects  and  re¬ 
fracts  the  light  rays  according  to  the  laws  of  optics.  No  absorption 
is  assumed  in  the  atmosphere,  unless  otherwise  mentioned. 

The  downward  radiance  just  below  the  ocean  surface  as  a  function 
of  the  zenith  angle  of  observation  is  shown  in  Fig.  2  for  the  case  of 
a  vanishingly  small  atmosphere  (t  tm 


=  0.0001)  and  a  solar  zenith  angle. 
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e0  =  11.44°.  Calculations  done  for  an  atmosphere  with  a  hundred-fold 
smaller  thickness  showed  no  differences  on  the  scale  of  these  figures. 
The  four  curves  are  for  ocean  albedo  values  of  w  *  1 ,  0.9,  0.5,  0.1. 

The  diffuse  radiance  shown  here  is  for  photons  that  have  undergone  one 
or  more  scattering  events;  thus  the  direct  solar  beam  is  not  shown.  In 
all  cases  the  radiance  increases  from  the  zenith  toward  the  critical 
angle  and  then  exhibits  only  a  slight  variation  from  the  critical  angle 
to  the  horizon.  The  upwelling  radiance  is  nearly  constant.  The  down- 
welling  radiance  necessarily  has  the  same  variation  with  angle  as  the 
upwelling  in  the  region  from  the  critical  angle  to  the  horizon,  since 
it  is  derived  by  total  internal  reflection  from  the  upwelling  in  this 
region.  The  downwelling  radiance  at  angles  less  than  the  critical 
angle  agrees  with  the  upwelling  radiance  at  the  same  angle  times  the 
reflectance  at  the  surface  appropriate  for  that  angle.  There  is  es¬ 
sentially  no  incoming  sky  radiation  from  this  very  thin  atmosphere. 

How  do  these  curves  change  as  the  atmosphere  becomes  thicker? 

The  curves  for  xatm  =  0.01  are  given  in  Fig.  3,  for  u)Q  =  1 ,  0.999,  0.99, 
0.9,  0.5,  0.1.  The  curves  are  essentially  unchanged  at  angles  greater 
than  the  critical  angle  from  Fig.  2.  However,  the  downward  radiance 
is  appreciably  greater  in  all  cases  for  angles  less  than  the  critical 
angle.  This  atmosphere  is  thick  enough  to  develop  some  sky  radiation 
which  enters  the  ocean  at  all  angles  less  than  the  critical  angle.  This 
sky  radiation  is  more  important  when  the  ocean  albedo  is  small,  since 
the  upwelling  radiance  in  the  ocean  is  smaller  in  this  case  also.  It 
is  interesting  that  a  spike  develops  in  the  radiance  at  the  critical 
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angle  when  u)o  =  0.1.  The  maximum  sky  radiance  develops  near  the  horizon 
in  this  thin  atmosphere  and  upon  entering  the  water,  it  Is  refracted 
into  a  narrow  range  of  angles  around  the  critical  angle. 

Curves  for  ratm  =  0.5  (corresponding  approximately  to  our  atmosphere 
in  the  visible)  are  given  in  Fig.  4.  The  upper  set  is  for  a  solar  zenith 
angle  eQ  =  11.44°,  while  the  lower  set  is  for  eQ  =  88.86°.  In  each  case 
curves  are  shown  for  o>0  =  1 ,  0.9,  0.7,  0.5,  0.1.  When  u>Q  <  0.8  (as  is 
the  case  for  the  actual  ocean),  these  curves  have  the  typical  shape  ob¬ 
served  with  the  largest  radiance  value  within  the  critical  cone.  When 
u>0  =  1  and  eQ  -  11.44°,  the  radiance  is  lower  within  the  critical  cone 
than  it  is  outside.  When  there  is  no  absorption  in  the  ocean,  the  up- 
welling  radiation  is  relatively  large  and  there  is  not  enough  sky  radia¬ 
tion,  even  with  this  atmospheric  thickness,  to  make  the  downwelling 
radiance  larger  within  the  critical  cone  than  it  is  outside. 

It  is  interesting  to  compare  the  lower  set  of  five  curves  for 
eQ  =  88.86°  with  the  upper  set  for  eQ  =  11.44°.  In  all  cases  there 
is  a  greater  drop  in  the  radiance  value  as  e  passes  through  the  critical 
angle  when  the  sun  is  near  the  horizon  than  when  it  is  near  the  zenith. 
Even  when  u)o  =  1  there  is  a  maximum  in  the  radiance  curve  just  before 
the  critical  angle  when  the  sun  is  near  the  horizon.  These  changes 
occur  because  there  is  relatively  more  sky  radiation  compared  to  the 
upwelling  radiation  in  the  ocean  when  the  sun  is  near  the  horizon. 

Similar  curves  are  given  in  Fig.  5  for  a  relatively  thick  atmosphere, 
Tatm  =  ^he  downwelling  radiance  just  below  the  ocean  surface  is 
nearly  constant  when  u>0  =  1 .  The  atmosphere  is  so  thick  that  the 


radiance  is  nearly  constant  above  the  ocean  surface.  Since  there  Is  no 
absorption  within  the  ocean,  this  nearly  constant  radiance  distribution 
holds  down  to  great  depths  until  it  is  eventually  influenced  by  the  nature 
of  the  bottom  surface  of  the  ocean. 

For  other  values  of  the  ocean  albedo  the  downwelling  radiance  just 
below  the  ocean  surface  has  a  maximum  at  the  zenith  and  decreases  ap¬ 
preciably  in  value  near  the  critical  angle. 

The  model  used  for  all  the  previous  results  has  no  absorption  in 
the  atmosphere.  It  is  an  interesting  theoretical  question  to  ask  how 
a  water  surface  changes  the  radiance  distribution  in  an  otherwise  homo¬ 
geneous  medium.  The  curves  in  Fig.  6  were  calculated  assuming  that 
u>0(atm)  =  w0(ocean),  i.e.  the  medium  above  and  below  the  interface  has 
the  same  properties.  The  optical  depth  of  the  interface  is  at  =  0.5, 
and  eQ  =  11.44°  and  88.86°.  A  comparison  of  this  figure  with  Fig.  4 
shows  a  much  larger  variation  in  the  radiance  near  the  critical  angle 
in  Fig.  4  than  in  Fig.  6  for  a  case  such  as  =  0.1  and  e  =  11.44°. 

This  is  caused  by  the  much  smaller  contribution  of  the  sky  radiation 
when  wo(atm)  =0.1  than  when  it  is  unity.  The  shape  of  the  curves  is 
much  more  nearly  similar  when  Figs.  4  and  6  are  compared  for  the  cases 
with  eQ  =  88.86°.  We  shall  return  to  the  question  of  the  influence  of  the 
interface  in  the  next  section. 

The  results  of  this  section  show  that  the  observed  decrease  in 
the  downwelling  radiance  just  below  the  ocean  surface  as  the  angle  of 
observation  Increases  through  the  critical  angle  only  occurs  when  there 
is  an  appropriate  balance  between  several  factors:  1.  the  optical 


thickness  of  the  atmosphere  must  be  large  enough  to  create  multiple 
scattered  sky  radiation  at  all  viewing  angles;  2.  the  single  scattering 
albedo  of  the  water  must  be  relatively  small  (wQ  <  0.7  approximately) 
so  that  the  upwelllng  radiance  just  below  the  ocean  surface  Is  a  small 
fraction  of  the  downwelllng  radiance  entering  the  ocean. 

V.  Dependence  of  Radiance  on  Depth  in  Medium 

The  interface  between  the  atmosphere  and  ocean  causes  Important 
changes  in  the  variation  of  radiance  with  zenith  angle,  especially  near 
the  interface.  In  this  section  various  examples  are  presented  which 
illustrate  how  the  radiance  distribution  varies  with  depth. 

In  order  to  study  the  change  in  shape  of  the  radiance  curves  with 
optical  depth,  it  is  convenient  to  plot  the  diffuse  downward  radiance 
at  a  particular  value  of  t  divided  by  the  diffuse  downward  irradiance 
for  the  same  value  of  t.  When  this  is  done,  it  is  possible  to  observe 
conveniently  the  change  in  shape  of  the  curve  on  one  graph;  otherwise 
the  simple  radiance  curves  are  eventually  dominated  by  the  exponential 
decrease  of  the  radiation  and  are  difficult  to  plot  on  a  single  graph. 

In  all  figures  in  this  section  the  radiance  has  been  divided  by  the 
diffuse  irradiance  at  that  optical  depth.  In  order  to  save  words  we 
refer  to  plots  of  the  "radiance  ratio"  understanding  that  the  ratio 
defined  above  is  the  actual  quantity  plotted. 

The  downward  radiance  ratio  is  shown  In  Fig.  7  when  the  sun  is  near 
the  horizon,  eo  =  88.86°,  the  ocean-atmosphere  Interface  is  at  an  optical 
depth  t,j  =  0.5,  u)Q(atm)  =  1,  wo(ocean)  =  0.5.  Near  the  top  of  the  at¬ 
mosphere  the  radiance  increases  toward  the  horizon  and  may  reach  a  maximum 
value  near  the  horizon  as  seen  in  the  curves  for  i*  0.1  and  0.5  (just 
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above  interface).  The  s^lld  curve  shows  the  downward  radiance  just  below 
the  Interface  with  the  typical  shape  increasing  from  the  zenith  to  near 
the  critical  angle  and  then  decreasing  by  an  order  of  magnitude  with  a 
slow  Increase  to  the  horizon.  The  sudden  decrease  near  the  critical 
angle  gradually  disappears  with  depth,  but  is  still  quite  evident  at 
t  *  1.5  (as  measured  from  the  top  of  the  atmosphere).  This  feature  has 
disappeared  when  t  =  10.5.  At  a  depth  of  t  *  50.5  the  radiance  ratio 
has  approached  the  asymptotic  radiance  distribution  closely,  at  least 
as  far  as  can  be  seen  on  the  scale  of  this  figure.  This  radiance  ratio 
distribution  at  great  depths  is  identical  to  the  one  obtained  for  a  homo¬ 
geneous  medium,  as  the  interface  no  longer  has  any  Influence  at  such 
depths. 

Only  one  parameter  is  changed  In  Fig.  8;  the  single  scattering  albedo 

In  the  ocean,  WQ(ocean)  s  0.1.  The  radiance  ratio  is  nearly  the  same  in 

the  atmosphere  at  optical  depths  down  to  the  Interface.  The  radiance 

ratio  just  below  the  interface  (solid  curve)  has  a  sharp  drop  of  over 

two  orders  of  magnitude  as  it  passes  through  the  critical  angle.  This 

drop  Is  still  quite  evident  at  t  *  1.5,  but  Is  beginning  to  disappear 

at  t  =  5.5.  At  still  larger  optical  depths  the  radiance  ratio  approaches 

its  asymptotic  form;  it  is  interesting  to  note  that  there  are  still 

noti cable  differences  between  the  actual  radiance  distribution  at  x  * 

100.5  and  the  asymptotic  distribution.  When  u>0  *  0.1,  extremely  large 

8 

optical  depths  are  required  to  establish  the  asymptotic  distribution. 

The  remaining  figures  in  this  section.  Figs.  9  -  12,  are  for  the 
sun  near  the  zenith,  eo  a  11.44°.  Each  figure  has  an  upper  and  lower 
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set  of  curves.  The  upper  set  is  for  a  homogeneous  medium  with  the 
indicated  single  scattering  albedo.  The  lower  set  is  for  a  medium 
with  an  atmosphere-ocean  interface  with  the  appropriate  reflection  and 
refraction  properties  for  an  upper  medium  with  an  index  of  refraction, 
n  =  1,  and  the  lower  medium  with  n  =  1.338.  We  are  Interested  in  studying 
the  effect  on  the  radiance  distribution  of  the  introduction  of  such  an 
Interface  into  a  homogeneous  medium;  thus  in  each  case  the  single  scat¬ 
tering  albedo  in  the  upper  medium  is  taken  as  equal  to  that  in  the  lower. 
In  all  these  figures  the  interface  is  taken  at  x^  =  0.5. 

In  Fig.  9,  u>Q  =  0.99.  The  downward  radiance  ratio  has  a  maximum 
at  the  horizon  for  small  optical  depths.  This  changes  into  a  distribution 
with  a  maximum  at  the  zenith  as  x  increases.  The  asymptotic  distri¬ 
bution  is  reached  on  the  scale  of  this  figure  at  x  ■  5.  The  radiance 
ratio  at  intermediate  optical  depths  is  quite  different  for  the  two 
cases  with  and  without  an  interface.  When  an  interface  is  present,  the 
radiance  ratio  just  below  the  interface  rises  toward  the  critical  angle 
and  then  shows  only  a  slight  increase  from  the  critical  angle  to  the 
horizon.  This  increase  still  occurs  at  x  =  1.5,  but  the  asymptotic  dis¬ 
tribution  is  reached  at  x  =  5.5,  the  same  distribution  as  for  the  homo¬ 
geneous  medium  without  interface. 

The  curves  in  Figs.  10  and  11  are  for  the  same  parameters  except 
that  u>0  =  0.5  and  0.1  respectively.  The  radiance  ratio  at  x  =  50  is 
still  somewhat  different  from  the  asymptotic  distribution  for  «  =  0.5. 
When  u>Q  =  0.1,  the  radiance  ratio  at  x  *  50  still  differs  from  the 
asymptotic  distribution  by  factors  of  more  than  three  at  many  angles. 


This  illustrates  the  general  fact  that  the  asymptotic  distribution  Is 

O 

valid  only  at  increasing  optical  depths  as  «0  becomes  smaller  . 

The  radiance  just  below  the  interface  for  both  w0  =  0.5  and  0.1 
increases  from  the  zenith  toward  the  critical  angle,  then  decreases 
rapidly  through  the  critical  angle,  and  finally  increases  slowly  toward 
the  horizon.  When  t  ^2.5,  the  radiance  ratio  is  a  monotonically  de¬ 
creasing  function  of  zenith  angle. 

In  general  the  shape  of  the  downward  radiance  ratio  curves  from  the 
top  of  the  atmosphere  to  the  top  of  the  interface  is  only  slightly  in¬ 
fluenced  by  whether  an  interface  is  present.  The  downward  radiance 
ratio  at  optical  depths  just  below  the  interface  is  naturally  influenced 
by  whether  or  not  an  interface  is  present.  This  influence  decreases  as 
t  increases  and  has  nearly  vanished  in  all  cases  when  t  =  5. 

A  similar  set  of  curves  for  the  upward  radiance  ratio  is  given  in 
Fig.  12  when  u)Q  =  0.5.  With  the  interface  at  =  0.5,  it  has  little 
influence  on  the  shape  of  the  upward  radiance  curves  when  t  >  0.5. 

However  the  curves  for  the  two  cases  with  and  without  the  interface  are 
quite  different  at  optical  depths  such  that  t  <  0.5.  At  the  top  of  the 
atmosphere  the  upward  radiance  ratio  is  much  stronger  at  the  horizon  and 
considerably  weaker  at  the  nadir  than  when  the  interface  is  present.  This 
is  due  to  two  factors:  the  stronger  reflection  of  the  direct  beam  from 
the  interface  for  angles  near  the  horizon  than  for  those  near  the  nadir 
and  the  fact  that  the  interface  prevents  an  appreciable  fraction  of  the 
upwelling  light  from  optical  depths  t  >  0.5  from  escaping  through  the 
interface. 


VI.  Dependence  of  Radiance  on  Position  of  Interface 

The  variation  of  the  radiance  with  the  position  of  the  atmosphere- 
ocean  Interface  is  Illustrated  by  the  Figs.  13-18.  The  radiance  just 
below  and  just  above  the  interface  is  given  in  Figs.  13-14  when  <»>0(atm)  = 

1  and  Wo(ocean)  =  0.5.  The  solar  zenith  angle  is  11.44°.  When  the  optical 
depth  of  the  interface  is  small  (x^  =  10’®),  the  downward  radiance  in¬ 
creases  from  a  very  small  (0.0054)  value  at  the  zenith  to  a  value  almost 
two  orders  of  magnitude  larger  near  the  critical  angle;  It  is  then  much 
more  nearly  constant  to  the  horizon.  The  upward  radiance  necessarily 
has  the  same  value  as  the  downward  at  angles  greater  than  the  critical 
angle  (because  when  it  is  totally  internally  reflected  it  becomes  the 
downward  radiance  at  these  angles);  the  upward  radiance  only  shows  minor 
variations  from  the  nadir  to  the  critical  angle.  As  x^  increases,  there 
is  the  same  qualitative  variation  of  the  radiance  for  x^  =  0.01,  but  when 
xj  ■  0.5,  the  incoming  radiation  from  the  sky  has  become  large  enough  to 
make  the  downward  radiance  just  below  the  surface  larger  for  angles  within 
the  critical  angle  than  without. 

These  results  may  be  compared  with  the  downward  and  upward  radiance 

just  above  the  surface  as  shown  in  Fig.  14.  The  downward  radiance  (sky 

radiation)  just  above  the  interface  Increases  uniformly  from  the  zenith 

-6  -2 

to  the  horizon  for  x^  =  10  and  10  .  When  x..  =  0.5  and  1  It  has  a 

maximum  near  the  horizon,  while  it  decreases  from  the  zenith  to  the 
horizon  when  x^  =  5.  On  the  other  hand,  the  upward  radiance  just  above 
the  interface  increases  from  the  nadir  to  the  horizon  for  x^  *  10*®  and 
10’^,  but  when  x^  =  0.5  it  passes  through  a  minimum  near  83°  and  then 
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increases  to  the  horizon.  This  increase  near  the  horizon  is  caused  by 
the  strong  reflection  of  the  downward  radiance  by  the  interface  for  in¬ 
cident  angles  near  90°.  When  x,.  =  1  and  5,  the  upward  radiance  just 
above  the  interface  again  increases  uniformly  from  the  nadir  to  the 
horizon.  Of  course  the  upward  and  downward  radiances  must  have  the 
same  value  at  the  horizon. 

In  order  to  determine  how  the  radiance  varies  with  solar  angle, 

the  downward  and  upward  radiance  just  below  the  surface  are  shown  in 

Fig.  15  for  the  same  parameters  as  in  the  last  two  figures,  except  that 

0Q  =  88.86°.  It  is  interesting  to  observe  that,  when  x^  has  the  very 
-4 

small  value  10  ,  the  downward  radiance  has  a  sharp  maximum  near  the 

_2 

critical  angle;  when  x.  =  10  ,  it  increases  to  a  maximum  near  the 
critical  angle  that  is  about  eight  times  greater  than  its  value  at 
the  zenith.  This  variation  develops  at  much  smaller  optical  depths 
when  the  sun  is  near  the  horizon  than  when  it  is  near  the  zenith  due 
to  the  much  greater  thickness  traversed  by  the  direct  solar  beam  in  the 
former  case. 

The  somewhat  complicated  interrelationships  between  the  different 
components  of  the  radiance  just  above  and  just  below  the  interface  are 
illustrated  in  Figs.  16,  17  and  18  for  models  with  u>o(atm)  =  <*>0(ocean)  = 
0.5,  0.9  and  0.1  respectively.  The  interface  is  taken  at  x..  =  0.5.  In 
these  figures  the  upward  and  downward  radiance  is  shown  both  just  above 
and  below  the  interface.  In  addition,  the  downward  radiance  just  above 
the  interface  is  shown  multiplied  by  the  reflection  between  air  and  water, 
R(A,  W),  (this  component  becomes  part  of  the  upward  radiance  just  above 
the  interface),  and  multiplied  by  the  transmission  of  the  interface  going 
from  air  to  water,  T  (A,  W)  (this  becomes  part  of  the  downward  radiance 
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just  below  the  interface).  In  the  latter  case  it  is  plotted  as  this 
part  of  the  radiance  would  appear  in  the  water  after  passing  through 
the  interface.  This  shows  that  the  reflected  downward  radiance  makes 
a  major  contribution  to  the  upward  radiance  just  above  the  interface 
only  at  angles  near  the  horizon,  because  the  reflectivity  of  water  is 
near  unity  only  at  such  angles.  Similarly  the  downward  radiance  trans¬ 
mitted  through  the  interface  is  the  major  component  of  the  downward 
radiance  just  below  the  interface  for  angles  from  the  zenith  to  those 
quite  near  the  critical  angle.  Of  course  it  can  make  no  contribution 
at  angles  greater  than  the  critical  angle. 

Also  in  this  figure  the  upward  radiance  just  below  the  interface 
is  shown  multiplied  by  the  reflection  between  water  and  air,  R(W,  A) 

(this  component  becomes  part  of  the  downward  radiance  just  below  the 
interface),  and  multiplied  by  the  transmission  of  the  interface  going 
from  water  to  air,  T(W,  A)  (this  becomes  part  of  the  upward  radiance 
just  above  the  interface).  This  is  plotted  as  this  part  of  the  radiance 
would  appear  in  the  atmosphere  after  passing  through  the  interface.  The 
reflected  upward  radiance  just  below  the  interface  becomes  the  entire 
downward  radiance  just  below  the  interface  at  angles  greater  than  the 
critical  angle,  but  is  a  major  part  of  the  downward  radiance  only  over 
a  small  additional  range  of  angles  slightly  smaller  than  the  critical 
angle.  The  transmitted  upward  radiance  that  passes  through  the  interface 
is  the  major  component  of  the  upward  radiance  above  the  interface  from 
the  nadir  to  angles  relatively  near  the  horizon. 

The  large  variation  in  the  shape  of  the  curve  representing  the 
downward  radiance  just  below  the  interface  as  uQ  and  -t.  vary  is  illu- 


strated  in  Fig.  19.  The  ratio  of  the  downward  radiance  at  8.52°  to  that 
at  49.37°  is  shown  as  a  function  of  w0.  The  first  angle  chosen  is  the 
quadrature  point  nearest  the  zenith,  while  the  second  is  the  first  quad¬ 
rature  point  past  the  critical  angle.  Thus  this  ratio  indicates  whether 
the  radiance  near  the  zenith  is  greater  than  or  less  than  the  radiance 
just  beyond  the  critical  angle.  When  eQ  =  11.46°  and  t.  =  0.0001,  this 
ratio  is  of  the  order  of  0.02,  but  as  the  interface  is  moved  down,  the 
ratio  rapidly  increases.  When  x.  =  5,  this  ratio  is  greater  than  100 
when  (joq  =  0.1,  indicating  that  the  radiance  is  much  larger  near  the 
zenith  than  at  the  horizon,  and  that  it  drops  appreciably  near  the 
critical  angle. 

In  Fig.  20  the  ratio  plotted  is  that  of  the  maximum  value  of  the 

downward  radiance  just  below  the  interface  to  the  value  at  49.37°.  These 

curves  indicate  that  this  ratio  increases  as  x.  and  e  increase  and  u>_ 

10  o 

decreases. 

VII.  Vector  and  Scalar  Irradiance 

The  vector  irradiance  or  flux,  H,  is  defined  as  the  integral  of 
the  product  of  the  radiance  with  the  cosine  of  the  angle  of  observation. 

The  upward  vector  irradiance,  Hu,  is  obtained  by  performing  the  integral 
over  all  angles  for  the  upward  hemisphere;  the  downward  vector  irradiance, 
H^,  is  obtained  from  the  integral  over  the  downward  hemisphere.  The  total 
vector  irradiance  is  the  integral  over  all  solid  angles  and  equals  Hy  -  H^. 
The  scalar  irradiance,  h,  is  the  Integral  of  the  radiance  over  solid 
angle  (upward  hemisphere,  downward  hemisphere,  or  all  solid  angles). 


The  reflectance  at  a  given  depth  in  the  medium  is  defined  as  R  =  Hu/H 
while  the  distribution  function  D  =  h/H  (defined  for  either  downward  or 
upward  irradiances) .  These  quantities  provide  a  relatively  simple  descrip 
tion  of  the  variation  of  the  radiation  with  depth  in  the  medium  without 
the  additional  complexity  of  having  to  specify  in  detail  the  variation 
of  the  radiance  with  angle. 

In  Table  I  the  total  upward  irradiance  including  the  direct  solar 

beam  reflected  from  the  water  interface  is  tabulated  in  the  second  and 

third  columns  as  a  function  of  the  optical  depth  from  the  top  of  the 

medium.  The  values  in  these  tables  are  for  e  =  11.44°  and  m  =  0.5  or 

o  o 

0.1.  The  values  in  the  columns  marked  "homog."  are  for  a  homogeneous 
medium  with  these  properties  and  no  interface.  When  an  interface  is 
present,  u>0(atm)  =  wQ(ocean)  =  0.5  or  0.1.  The  values  in  the  columns 
marked  "with  interface"  have  an  atmosphere-water  interface  at  an  optical 
depth  of  0.5.  Thus  a  comparison  of  these  two  columns  shows  directly 
the  influence  of  the  interface. 

When  u)o  =  0.5  the  upward  irradiance  is  smaller  above  the  interface 

( x  <  0.5)  in  the  medium  with  the  interface  than  at  comparable  depths  in 

the  homogeneous  medium,  while  it  is  greater  at  corresponding  depths  in 

the  medium  with  the  interface  when  t>0.5.  This  variation  is  in  the 

opposite  sense  when  w  =0.1.  In  the  latter  case  there  is  so  much  ab- 

o 

sorption  in  the  medium  that  the  radiation  from  the  direct  beam  reflected 
into  the  upward  direction  at  the  interface  makes  an  appreciable  contribu¬ 
tion  to  the  irradiance  when  t  <0.5.  On  the  other  hand  when  there  Is  less 
absorption,  toQ  =  0.5,  the  upward  irradiance  is  much  larger  at  all  levels 
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compared  to  the  case  with  uQ  =  0.1.  Only  part  of  the  upward  irradiance 
just  below  the  interface  is  transmitted  through  the  interface,  with 
the  result  that  the  total  upward  irradiance  is  less  when  x<0.5  for  the 
medium  with  the  interface  than  for  the  homogeneous  medium.  The  reflected 
solar  beam  makes  only  a  minor  contribution  to  the  total  upward  irradiance 
above  the  interface  when  iuo  =  0.5. 

The  fourth  and  fifth  columns  in  Table  I  give  the  total  downward  ir¬ 
radiance  including  the  direct  solar  beam  for  these  two  media.  In  this 
study  the  downward  flux  is  always  normalized  to  unity  for  a  plane  per¬ 
pendicular  to  the  incoming  beam.  When  ojq  =  0.5,  the  downward  irradiance 
is  greater  just  below  the  interface  than  it  is  in  the  homogeneous  medium 
at  the  same  depth  due  to  the  upwelling  radiation  reflected  into  the 
downward  direction  by  the  interface.  The  opposite  situation  prevails 
when  u>o  =  0.1,  since  the  upwelling  irradiance  is  so  small;  in  this  case 
the  loss  of  some  of  the  downward  irradiance  by  reflection  at  the  inter¬ 
face  is  the  determining  factor  in  this  variation. 

The  last  two  columns  give  the  total  scalar  irradiance  (hu  +  h^) 
which  is  proportional  to  the  heating  rate  at  the  given  depth.  This 
quantity  is  less  for  the  medium  with  the  interface  than  for  the  homo¬ 
geneous  medium  when  t<0.5  and  uiQ  =  0.5,  while  it  is  larger  when  t>0.5. 
There  is  the  opposite  variation  when  =  0.1.  All  of  the  quantities  in 
Table  I  decrease  nearly  exponentially  with  optical  depth  when  measured 
reasonably  far  from  boundaries. 

On  the  other  hand  the  quantities  in  Table  II  are  remarkably  constant 
throughout  the  medium.  The  second  and  third  columns  of  Table  II  give  the 
reflectance.  The  fourth  and  fifth  columns  give  the  ratio  of  the  total 
upward  irradiance  to  the  total  downward  irradiance  (total  means  that 


the  direct  beam  is  included  in  the  downward  irradiance  and  that  the  re¬ 
flected  direct  beam  from  the  interface  is  included  in  the  upward  irradiance 
diffuse  means  that  these  quantities  are  not  included). 

The  distribution  function  Du  =  hu/Hu  is  given  in  the  sixth  and  seventh 
columns.  In  this  expression  Hu  and  hu  are  the  total  upward  irradiance 
and  total  upward  scalar  irradiance  respectively,  including  the  contribu¬ 
tion  from  the  solar  beam  reflected  at  the  interface.  The  eighth  and 
ninth  columns  give  the  same  quantity  for  the  downward  radiation,  = 
h^/H^,  with  the  use  of  the  total  irradiance  and  scalar  irradiance,  in¬ 
cluding  the  contribution  from  the  solar  beam.  The  last  two  columns  give 
the  distribution  function  for  the  downward  radiation,  but  only  the  diffuse 
contribution  to  the  irradiance  and  scalar  irradiance  is  used. 

The  distribution  function  is  important  because  it  depends  on  the 
variation  of  the  radiance  with  angle.  If  the  radiation  is  monodirectional 
with  direction  cosine  then  D  =  l/wQ,  while  if  the  radiance  is  uniformly 
distributed  with  angle,  then  0  =  2.  If  0  is  close  to  unity,  the  radiance 
has  a  maximum  near  the  zenith  (for  downwelling  radiation)  and  only  small 
values  near  the  horizon.  Similarly  if  D  >  2,  the  radiance  is  greater  near 
the  horizon  than  near  the  zenith.  Thus  in  the  last  column  of  Table  II, 

D  =  2.943  in  the  medium  with  an  interface,  when  t  =  0.1  and  =  0.5,  since 
the  diffuse  radiation  develops  first  near  the  horizon  at  small  optical 
depths  because  of  the  greater  effective  optical  thickness  there.  Note 
the  different  values  of  D  just  above  and  just  below  the  interface.  Just 
below  the  interface,  D  =  1.863,  since  the  downwelling  radiation  is  largely 
confined  within  the  critical  cone  around  the  zenith  direction.  At  great 
depths  D  =  1.359  corresponding  to  the  radiance  distribution  peaked  around 
the  zenith  that  develops  in  the  diffusion  region  in  an  absorbing  medium. 


The  variation  of  some  of  these  quantities  is  shown  in  Fig.  21  for 
the  case  eQ  =  11.44°,  x^  =  0.5,  ai0(atm)  =  1,  ^(ocean)  =  0.5.  The 
following  quantities  are  indicated  as  a  function  of  the  optical  depth 
from  the  top  of  the  atmosphere:  total  and  diffuse  upward  and  irradiance, 
total  and  diffuse  downward  irradiance,  the  net  flux  equal  to  the  dif¬ 
ference  between  the  total  downward  and  total  upward  irradiances,  and 
the  derivative  with  respect  to  optical  depth  of  the  total  flux  (which 
is  proportional  to  the  heating  rate). 

Under  what  conditions  does  the  diffuse  flux  increase  or  decrease 
upon  passing  through  the  interface?  This  is  illustrated  in  Table  III 
which  gives  the  ratio  of  the  downward  diffuse  irradiance  (direct  solar 
beam  not  included)  just  below  the  interface  to  that  just  above;  the 
same  ratio  for  the  upward  diffuse  irradiance  (direct  solar  beam  reflected 
by  interface  not  included)  is  also  shown.  In  most  cases  given  in  the 
Table,  the  downward  diffuse  irradiance  is  greater  just  below  the  inter¬ 
face  than  above.  In  some  cases  it  is  two  to  three  or^rs  of  magnitude 
greater,  e.g.  when  the  interface  is  near  the  top  of  the  atmosphere  and 
the  sun  is  near  the  zenith.  In  these  cases  there  is  insufficient  at¬ 
mosphere  for  any  appreciable  downwelling  scattered  radiation  to  develop. 
The  downward  diffuse  radiance  just  below  the  interface  is  less  than  that 
just  above  only  in  cases  where  the  absorption  in  the  ocean  is  appreciable 
and  the  interface  is  at  an  intermediate  optical  depth.  When  the  ocean 
absorption  is  large,  there  is  a  much  smaller  upwelling  irradiance  in  the 
ocean  and  thus  relatively  little  radiation  is  reflected  at  the  interface 
into  the  downward  direction.  The  ratios  for  the  upward  irradiance  can 
be  explained  in  a  similar  manner. 
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VIII.  Effect  of  Phase  Function 

All  results  given  in  the  previous  sections  assume  a  Rayleigh  phase 
function  for  the  scattering  centers  in  both  the  atmosphere  and  ocean. 

It  is  realized,  of  course,  that  the  results  would  be  changed  for  other 
phase  functions.  In  order  to  explore  how  large  the  changes  are  that 
might  be  anticipated,  some  of  the  calculations  were  repeated  with  a  more 


realistic  phase  function  in  the  ocean.  The  phase  function  chosen  was 


proportional  to  A(w  +  0.75)"^  +  B,  where  A  =  2.52  x  10"^,  B  =  0.1,  and 


u  is  the  cosine  of  the  scattering  angle.  The  ratio  of  scattering  into 
the  forward  hemisphere  to  that  into  the  backward  is  9,  while  the  ratio 


of  the  scattering  at  0°  to  that  at  180°  is  152.  The  form  of  this  function 


was  chosen  merely  because  it  could  be  included  in  the  computer  program 
with  relative  ease.  However,  it  does  reproduce  some  of  the  features  of 
a  realistic  ocean  phase  function,  such  as  a  strong  forward  to  backward 
scattering  ratio  and  a  peak  around  0°. 

The  downward  and  upward  radiance  just  below  the  interface  is  shown 
in  Fig.  22.  In  order  to  study  the  transition  from  Rayleigh  to  a  Mie-type 
phase  function,  calculations  were  made  for  several  values  of  the  parameter 
6  =  °5R/as7>  where  o$R  is  the  scattering  cross  section  in  the  ocean  for 
Rayleigh  scattering  processes  and  is  the  total  scattering  cross  section 
from  all  processes.  Thus  when  6=1,  the  scattering  is  entirely  Rayleigh 
and  the  results  are  the  same  as  those  given  in  the  previous  sections.  When 
6=0,  there  is  not  Rayleigh  scattering  and  all  scattering  is  by  the  phase 
function  given  in  the  previous  paragraph. 


The  upwelling  radiance  shown  in  Fig.  22  shows  that  this  quantity  de¬ 
creases  appreciably  as  the  scattering  becomes  predominately  in  the  forward 
direction,  since  the  downwelling  photons  have  less  chance  to  be  scattered 


into  an  upward  direction.  This  has  relatively  little  effect  on  the 
downwelling  radiance  within  the  critical  cone.  However  the  effect  is 
large  outside  the  critical  cone,  since  the  downwelling  radiance  in  this 
range  of  angles  is  entirely  derived  from  the  upwelling  radiance  totally 
internally  reflected  at  the  atmosphere -water  interface. 

IX.  Conclusion 

The  influence  of  the  various  properties  of  an  atmosphere-ocean  system 
on  the  radiance  distribution  has  been  explored.  The  downwelling  radiance 
just  below  the  ocean  surface  has  a  well-known  maximum  within  the  cone 
extending  out  to  the  critical  angle  and  then  decreases  by  one  to  two 
orders  of  magnitude  toward  the  horizon.  This  observed  variation  in  the 
visible  in  actual  bodies  of  water  depends  not  only  on  the  laws  of  re¬ 
flection  and  refraction  at  the  interface,  but  crucially  on  two  other 
properties  of  the  medium.  The  atmosphere  must  have  sufficient  optical 
depth  to  develop  enough  multiple  scattered  sky  radiation  at  all  angles 
to  provide  the  source  for  the  enhanced  light  observed  within  the  critical 
angle  just  beneath  the  ocean  surface.  In  addition  the  water  must  have 
a  fairly  high  absorption  for  visible  light,  so  that  the  ratio  of  the 
upwelling  to  downwelling  light  within  the  ocean  is  small.  When  this 
happens  the  downwelling  light  outside  the  critical  angle  just  beneath 
the  ocean  surface  (derived  by  total  internal  reflection  from  the  weak 
upwelling  light)  is  one  or  two  orders  of  magnitude  less  than  the  down¬ 
welling  light  within  the  critical  angle.  Numerous  examples  given  here 
show  that  just  the  opposite  situation  occurs  when  these  conditions  are 
not  satisfied. 


1 
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The  radiance  distribution  with  depth  for  a  homogeneous  medium  is 
compared  with  that  for  an  identical  medium,  but  with  an  atmosphere -ocean 
interface.  The  influence  of  the  interface  on  the  radiance,  irradiance, 
heating  rate  and  many  other  properties  of  the  radiation  has  been  illu¬ 
strated.  The  variation  of  these  quantities  with  depth  is  complex  in 
its  dependence  on  the  scattering  and  absorption  properties  of  the  at¬ 
mosphere  and  ocean,  including  the  phase  function  for  scattering  as  well 
as  the  optical  depth  from  the  top  of  the  atmosphere  of  the  atmosphere- 
ocean  interface. 
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TABLE  I 


<*>o  =  0.5  ;  eQ  =  11.44° 

Hu(total)  Hd(total)  hu  +  hd 


T 

homog. 

with 

interface 

homog. 

with 

interface 

homog. 

with 

interface 

0 

0.1192 

0.1094 

0.9801 

.9801 

1.240 

1.220 

.1 

.1118 

.1004 

.9113 

.9108 

1.207 

1.182 

.5* 

8.346-2*** 

6.255-2 

.6673 

.6647 

0.9442 

0.8894 

.5** 

8.346-2 

8.598-2 

.6673 

.6882 

.9442 

.9682 

1 

5.604-2 

5.784-2 

.4422 

.4565 

.6494 

.6701 

2 

2.411-2 

2.480-2 

.1876 

.1932 

.2859 

.2935 

5 

1.665-3 

1.717-3 

1.274-2 

1.317-2 

2.021-2 

2.080-2 

10 

1.665-5 

1.734-5 

1  .256-4 

1.318-4 

2.033-4 

2.127-4 

20 

1.395-9 

1.494-9 

1.054-8 

1 .129-8 

1.724-8 

1.844-8 

50 

6.638-22 

7.238-22 

5.006-21 

5.458-21 

8.224-21 

8.966-21 

O 

1! 

o 

3 

;  eQ  =  H 

.44° 

0 

1.754-2 

2.304-2 

0.9801 

0.9801 

1 .035 

1.040 

.1 

1.594-2 

2.191-2 

.8893 

.8893 

0.9457 

0.9523 

.5* 

1.080-2 

1.886-2 

.6003 

.6004 

.6467 

.6524 

.5** 

1.080-2 

1.061-2 

.6003 

.5921 

.6467 

.6317 

1 

6.591-3 

6.510-3 

.3657 

.3623 

.3967 

.3898 

2 

2.436-3 

2.424-3 

.1349 

.1346 

.1472 

.1457 

5 

1.200-4 

1.220-4 

6.627-3 

6.755-3 

7.282-3 

7.368-3 

10 

7.738-7 

8.143-7 

4.269-5 

4.503-5 

4.704-5 

4.931-5 

20 

3.153-11 

3.511-11 

1.739-9 

1.939-9 

1.918-9 

2.129-9 

50 

4.776-26 

2.630-24 

1 .205-22 

1.452-22 

1.013-22 

1.597-22 

♦above  interface 
♦♦below  interface 
♦♦♦represents  8.346  x  10“* 
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TABLE  III 

Ratio  Diffuse  Irradiance  Just  Below  to  Diffuse  Irradiance  Just  Above 
Interface 


Downward  Upward 


Ti 

“o 

9  =11.44° 
0 

0  =76.28° 
0 

G  =  '8.86° 

0 

0  =11.44° 
0 

0  =76.28° 
0 

0  =88.86' 
0 

0.0001 

1 

6334 

2181 

24.89 

1 .896 

1.982 

1.937 

0.9 

3540 

900.7 

9.442 

1.968 

2.089 

2.048 

0.5 

1070 

236.7 

2.956 

1.924 

2.155 

1.968 

0.1 

155.5 

33.33 

1.101 

1.861 

2.146 

1.231 

0.01 

1 

64.28 

22.83 

1.796 

1.893 

1.939 

1.649 

0.9 

36.01 

9.856 

1.155 

1.960 

2.045 

1.312 

0.5 

11.49 

3.185 

0.868 

1.903 

1.983 

0.611 

0.1 

2.397 

1 .147 

0.830 

1.730 

1.321 

0.111 

0.5 

1 

2.757 

1 .893 

1.770 

1.829 

1.778 

1.771 

0.9 

1.784 

1 .261 

1.210 

1 .814 

1.620 

1 .600 

0.7 

1.351 

1 .074 

1.056 

1.669 

1.315 

1.286 

0.5 

1.155 

0.999 

0.995 

1.472 

0.992 

0.964 

0.1 

0.952 

0.926 

0.937 

0.593 

0.231 

0.222 

1 

1 

2.156 

1.794 

1.784 

1.810 

1.782 

1.786 

0.9 

1.450 

1.223 

1 .221 

1.757 

1.636 

1.648 

0.5 

1.066 

1.004 

1.006 

1.316 

1.026 

1.050 

0.1 

0.950 

0.943 

0.948 

0.434 

0.248 

0.260 

5 

1 

1.794 

1.788 

1.788 

1.792 

1.791 

1.790 

0.9 

1 .230 

1.226 

1 .226 

1.674 

1.671 

1.672 

0.5 

1.013 

1.012 

1.012 

1.108 

1.101 

1.101 

0.1 

1 .000 

1.000 

1 .000 

0.290 

0.286 

0.286 
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Fig.  1. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


Fig.  5. 


Fig.  6. 


Fig.  7. 


Legends  for  Figures 


Downward  radiance  just  below  ocean  surface  as  function  of  zenith 

angle  of  observation  e  for  solar  zenith  angle  0Q  =  11.44°.  The 

top  curve  has  single  scattering  albedo  in  the  ocean,  u>Q  =  0.5 

and  an  atmosphere  of  optical  thickness  xatm  =  0.5.  The  middle 

curve  has  &>_  =  0.5  and  t  .  =  0.0001.  The  bottom  curve  has 

o  atm 

“o  K0  and  Tatm  =  °*5* 


Downward  radiance  just  below  ocean  surface  as  function  of  e  for 
eQ  =  11.44°,  Tatm  =  0.0001  and  «  =  0.1,  0.5,  0.9,  1.0. 

Downward  radiance  just  below  ocean  surface  as  function  of  0  for 
0O  =  11.44°,  Tatm  =  0.01  and  u>0  =  0.1,  0.5,  0.9,  0.99,  0.999,  1.0. 

Downward  radiance  just  below  ocean  surface  as  function  of  0  for 
en  =  11.44°  and  88.86°,  x  .  =  0.5  and  a.  =  0.1,  0.5,  0.7,  0.9, 

o  a  cm  0 

1.0. 


Downward  radiance  just  below  ocean  surface  as  function  of  0  for 
eQ  =  11.44°  and  88.86°,  Tatm  =  5  and  u>0  =  0.1,  0.5,  0.9,  1.0. 

Downward  radiance  just  below  ocean  surface  as  a  function  of  0 
for  0Q  =  11.44°  and  88.86°,  x..  =  0.5,  u)Q(atm)  =  u>Q(ocean)  and 
u>0  =  0.1,  0.5,  0.7,  0.9,  0.99,  1.0. 

Downward  radiance  divided  by  diffuse  downward  scalar  irradiance 

as  a  function  of  0  for  e  =  88.86°,  x.  =  0.5,  u>  =  0.5  at  various 

o  i  o 

optical  depths  measured  from  top  of  atmosphere. 


Fig.  8.  Same  as  Fig.  7,  except  uo  =  0.1. 

Fig.  9.  Downward  radiance  divided  by  diffuse  downward  irradiance  as  a 

function  of  e  for  eQ  =  11.44°  at  various  optical  depths.  The 

upper  curves  are  for  a  homogeneous  medium  with  u)0  =  0.99.  The 

lower  curves  are  for  a  medium  with  an  interface  at  an  optical 

depth  T-  =  0.5  and  with  w  (atm)  =  w  (ocean)  =  0.99. 

1  o  o 

Fig.  10.  Same  as  Fig.  9  except  u)q  =  0.5  for  upper  curves  and  ioo(atm)  = 
u>0(ocean)  =  0.5  for  lower  curves. 

Fig.  11.  Same  as  Fig.  9  except  =  0.1  for  upper  curves  and  wo(atm)  = 

u  (ocean)  =  0.1  for  lower  curves, 
o 

Fig.  12.  Upward  radiance  divided  by  upward  irradiance  as  a  function  of 

e  for  eQ  =  11.44°  at  various  optical  depths.  The  upper  curves 

are  for  a  homogeneous  medium  with  =  0.5.  The  lower  curves 

are  for  a  medium  with  an  interface  at  an  optical  depth  t..  =0.5 

and  with  w  (atm)  =  u>  (ocean)  =  0.5. 
o  o 

Fig.  13.  The  upward  and  downward  radiance  (divided  in  both  cases  by  down¬ 
ward  radiance)  just  below  the  interface  for  oo  =  11.44°,  w0(atm) 
1,  ioQ(ocean)  =  0.5,  with  the  interface  at  the  following  optical 
depths:  ^  =  10*6,  10‘2,  0.5,  1,  5. 
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Fig.  14.  The  upward  and  downward  radiance  just  above  the  interface  for 

0Q  =  11.44°,  u)Q(atm)  =  1,  u)Q(ocean)  =  0.5,  with  the  interface 

-6  -4  _? 

at  the  following  optical  depths:  x..  =  10  ,10  ,  10  ,  0.5,  1,  5. 

Fig.  15.  The  upward  and  downward  radiance  just  below  the  interface  for 

0„  =  88.86°,  w  (atm)  =  1,  w  (ocean)  =  0.5,  with  the  interface  at 
0  0  0 

-6  -4  -? 

the  following  optical  depths:  =  10  ,10  ,10  ,  0.5,  1,  5. 

Fig.  16.  The  upward  and  downward  radiance  just  above  and  just  below  the 

interface  for  0Q  =  11.44°,  u)Q(atm)  =  u)Q(ocean)  =  0.5,  and  = 

0.5.  These  radiances  are  also  shown  multiplied  by  the  reflection 
and  transmission  of  the  interface  going  from  air  to  water,  R(A,  W) 
and  T(A,  W)  respectively.  They  are  also  shown  multiplied  by  the 
reflection  and  transmission  of  the  interface  going  from  water  to 
air,  R(W,  A)  and  T(W,  A)  respectively. 

Fig.  17.  Same  as  Fig.  16  except  wQ(atm)  =  u>o(ocean)  =  0.9. 

Fig.  18.  Same  as  Fig.  16  except  wo(atm)  =  u)Q(ocean)  =  0.1. 

Fig.  19.  Ratio  of  downward  diffuse  radiance  just  below  interface  at  8.52° 

to  that  at  49.37°  (just  beyond  critical  angle)  as  a  function  of 

u>  for  various  values  of  x.  and  0  . 
o  10 

Fig.  20.  Ratio  of  maximum  value  of  downward  diffuse  radiance  just  below 

interface  to  radiance  value  at  49.37°  (just  beyond  critical  angle) 

as  a  function  of  w  for  various  values  of  x.  and  e  . 

0  10 
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Fig.  21.  The  total  and  diffuse  upward  irradiance,  Hu>  the  total  and 

diffuse  downward  irradiance,  H^,  the  total  net  flux,  -  Hu, 

and  the  derivative  of  the  total  flux  with  respect  to  optical 

depth  (proportional  to  heating  rate)  as  a  function  of  optical 

depth  in  the  medium  for  eQ  =  11.44°,  t..  =  0.5,  <oo(atm)  =  1, 

o>  (ocean)  =  0.5. 
o 

Fig.  22.  The  upward  and  downward  radiance  just  below  the  interface  for 
various  values  of  8,  the  ratio  of  the  Rayleigh  scattering  cross 
section  to  the  total  (Rayleigh  plus  Mie  —  see  text)  scattering 
cross  section. 
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